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: Abstract 

We study the Standard Model prediction for the mass difference between the two neutral D 
meson mass eigenstates, Am. We derive a dispersion relation based on heavy quark effective 
theory that relates Am to an integral of the width difference of heavy mesons, Ar, over varying 
values of the heavy meson mass. Modeling the m£)-dependence of certain D decay partial widths, 
we investigate the effects of SU (3) breaking from phase space on the mass difference. We find that 
Am may be comparable in magnitude to Ar in the Standard Model. 



I. INTRODUCTION 



The mixing and decay of K, B, and D mesons are sensitive probes of physics beyond the 
Standard Model. Among the many processes that one might study, flavor-changing neutral 
current D decays and D° — D° mixing provide unique information, because in the Standard 
Model (SM) they occur via loop diagrams involving intermediate down-type quarks. In par- 
ticular, because of severe CKM and GIM suppressions, the mixing of D mesons is expected 
to be quite slow, and thus the D system is one of the most intriguing probes of new physics 
in low energy experiments [1]. 

We begin by recalling the formalism for heavy meson mixing. Using standard notation, 
the expansion of the off-diagonal terms in the neutral D mass matrix to second order in 
perturbation theory is given by 

2 /12 2m D V 2m D V m D -E n + ie yj 

The first term represents the AC = 2 contributions that are local at the scale // ~ mo. It 
contributes only to M 12 , and is expected to be very small unless it receives large enhancement 
from new physics. The second term in Eq. (1) comes form double insertion of AC = 1 
operators in the SM Lagrangian and it contributes to both M 12 and Ti 2 . It is dominated 
by the SM contributions even in the presence of new physics. Two physical parameters that 
characterize the mixing are 

Am Ar 

* = —, y=2r, ( 2 ) 

where Am and Ar are the mass and width differences of the two neutral D meson mass 
eigenstates and V is their average width. Because of the GIM mechanism the mixing am- 
plitude is proportional to differences of terms suppressed by m 2 , sb /m 2 v , and so D° — D° 
mixing is very slow in the SM [2]. The contribution of the b quark is further suppressed 
by the small CKM elements \V ub V* b \ 2 /\V us V c *\ 2 = £>(l(r 6 ), and can be neglected. Thus, 
the D system essentially involves only the first two generations, and therefore CP violation 
is absent both in the mixing amplitude and in the dominant tree-level decay amplitudes, 
and will be neglected hereafter. Once the contribution of b quarks is neglected, the mixing 
vanishes in the flavor SU(3) limit, and it only arises at second order in SU(3) breaking if 
SU (3) breaking can be treated analytically [3] 

x , y ~ sin 2 9 C x [SU(3) breaking] 2 , (3) 



where 8c is the Cabibbo angle. Precise calculations of x and y in the SM are not possible 
at present, because the charm mass is neither heavy enough to justify inclusive calculations, 
nor is it light enough to allow a few exclusive channels to give a reliable estimate. 

According to Eq. (3), computing x and y in the SM requires a calculation of SU(3) 
violation in the decay rates. There are many sources of SU(3) violation, most of them 
involving nonperturbative physics in an essential way. In Ref. [3], SU(3) breaking arising 
from phase space differences was studied; computing them in two-, three-, and four-body D 
decays, it was found that y could naturally be at the level of one percent. This result can 
be traced to the fact that the SU (3) cancellation between the contributions of members of 
the same multiplet can be badly broken when decays to the heaviest members of a multiplet 
have small or vanishing phase space. This effect is manifestly not included in the OPE-based 
calculations of D° — D° mixing, which cannot address threshold effects. 

The purpose of the present paper is to address the following question: if the dominant 
577(3) breaking mechanism is indeed the one studied in Ref. [3], and it gives rise to y at the 
percent level, then can x naturally be comparably large? This is particularly relevant because 
the present experimental upper bounds on x and y are at the few times 10~ 2 level [4, 5] and 
are expected to significantly improve (for a review of the experimental situation, see Ref. [6]). 
To interpret the results from future measurements of x and y, and possibly establish the 
presence of new physics, we need to know the allowed range in the SM. In particular, since 
new physics can only contribute to x, an experimental observation of x ^> y would imply 
a large new physics contribution to D° — D° mixing. Although y is determined by SM 
processes, its value still affects the sensitivity to new physics [7]. 

In this paper we study the SM predictions for x/y due to SU{3) breaking from final 
state phase space differences. In Sec. II we derive a dispersion relation using Heavy Quark 
Effective Theory (HQET) that relates Am to Ar. To compute Am, we need a calculation 
of Ar for varying heavy meson mass, so we review its calculation from Ref. [3] in Sec. III. 
In Sec. IV, we calculate Am and present numerical results. We find that despite the fact 
that SU (3) breaking in phase space affects x in a different way than it affects y, the final 
estimates of x and y are comparable. We present our conclusions in Sec. V and discuss the 
implications of our findings for experimental searches for new physics in D° — D° mixing. 
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FIG. 1: The correlator in Eq. (6). The black boxes denote the weak Hamiltonian, the wavy lines 
show external momenta inserted, and the gray area represents hadronic intermediate states. 

II. DERIVATION OF THE DISPERSION RELATION 

We start by reviewing the relevant formalism for D° — D° mixing. Equation (1) im- 
plies that the mass eigenstates are linear combinations of the weak interaction eigenstates, 
^1,2) = P =t <? Since we neglect the effects of intermediate states containing a b 

quark, |-D 12 ) are also CP eigenstates, CP\D±) = ±\D±). Their mass and width differences 
are 

Am = m D+ - m D _ = 2M 12 , Ar = T D+ - T D _ = 2T 12 . (4) 
Neglecting the small contribution from the local AC = 2 operators, Eq. (1) gives 

Am = 1 F y (D°\n w \n}(n\n w \D ) + (D°\H w \n)(n\H w \D } 
2m D ^ m D - E n 

Ar = J- \(D°\H w \n)(n\H w \D } + (D°\H w \n) (n\H w \D )] (2n)6(m D - E n ) , (5) 

where P denotes the principal value prescription, the sum is over all intermediate states, n, 
and it implicitly includes (2ti) 3 5 3 (pd —Pn)- 

To derive a dispersion relation between Am and AT, consider the following correlator 

S PD (g) = % Jd A z (D(p D )\T \H w {z) 7^(0)] \D(p D )) e^™> z . (6) 

Here po is a label given by the momentum of the on-shell D meson state (satisfying p 2 D = m 2 D ) 
and q—pD is an auxiliary four- vector that inserts external momentum to the weak interaction 
(see Fig. 1). There is no simple physical interpretation of £ except at q — p^, where H PD (p£,) 
is related to physical properties of D mesons. Inserting a complete set of states in Eq. (6) 
and comparing with Eq. (5), we find 

' T, PD {p D )= (Am-'- AT). (7) 



2m D PD ^ LJ ' V 2 
The correlator E PD (g) is an analytic function of q (but not of po) with a cut in the complex 



q° plane for q° > J \q \ 2 + 4m^ for a fixed q. 



To write the dispersion relation in terms of physical quantities, i.e., to give Ti PD (q) for 
q 7^ Pd a physical interpretation, we need to eliminate the heavy quark mass dependence 
from Eq. (6). 1 The momentum of a heavy meson H containing a heavy quark Q can be 
written as p^ = mnv^, with v 2 = 1. We can decompose Q as 



Q( z ) = e - tm ^ z hi Q) (z) + e +im ® v - z h ( Q\z) + . . . , 



(8) 



where the HQET fields and respectively annihilate a heavy Q quark and create 
a heavy Q antiquark with four- velocity v. Here and in the rest of this section the ellipses 
denote terms suppressed by a relative factor of Aqcd/^c- The AC = 1 weak Hamiltonian 
contributing to neutral D meson mixing is 
4G F 



V2 



(9) 



where qi^ = d or s, and the four-quark operators, suppressing their Dirac structures, are of 
the form 



Oi ~ qiq2uc = e 



-im c v-z 



q^util 



: im c v-z 



q\q2Uti\, 



uh^ + ... . 



(10) 



In Eq. (9) TC W contains the light quark fields, the Wilson coefficients, and summation over 
operators. We also replace the QCD states \D) by HQET states \H(v)), 



\D(p = m D v)) = y/mE\H(v)) + ... . 



(11) 



The new states have a normalization that is independent of the heavy quark mass [9] . Then 
Eq. (6) yields 



S PD (q)=im D J d 4 z (H(v) | T { [ e ~~* H w h® (z) + e m ^ z HJi^ ( 



x 



H w h^(0) + H w h ( c\0)] } e ^- p °> z \H(v)) + .... 

The only nonzero contributions to this correlator involve a single h and h field each, 

S PD (?)=mD Jd*z (H(v)\ {e~ m ^ z T[n w h^\z),Hjh^(0)_ 

H w h^\z),H w h^(0)]}e^™> z \H(v)) + .... 



(12) 



T 



(13) 



The two terms in Eq. (13) behave differently in the HQET limit m c — > oo with q fixed. 
The term proportional to exp[i(q— p D — m c v) • z] oscillates infinitely rapidly and is integrated 



1 The method of using HQET to derive a dispersion relation in the heavy quark mass was developed first 
in Ref. [8], where it was used to study the inclusive nonleptonic heavy meson decay rate. 
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out at the heavy scale. It should be removed from the effective theory and replaced by a 
local H^ c=2 contribution that can be included as a matrix element of AC = 2 operators. 
Such contributions are estimated to give rise to x and y at or below the 10~ 3 level [10-12], 2 
and since we are interested in the question whether x could be near the percent level, we 
can neglect them. 

By contrast, the term proportional to exp[i(q — po + rn c v) • z] becomes independent of 
m c as m c ^ oo. Recalling that pr> = rriDV, we have 

S PD (g) = im D j d 4 z (H(v)\T[nJii c \z),H w hi c \0)] e^~ Kv > z \H(v)) + . . . , (14) 

where A = mp — m c + (9(AQ CD /m c ). It is convenient to define 

=i [d 4 z(H(v)\T\n w h^\z),n w h^\0)]e i ^^ (15) 



which is manifestly independent of the heavy quark mass. It follows that 

S pd(^) = m D i: v (q) + ... , 
and Eq. (7) becomes to leading order in Aqcd/^c 

= -2 Am(E) + i AT(E) , 



(16) 



(17) 



where E = ^fq 1 , and Am(E) and AT(E) can be interpreted as the mass and the width 
differences of neutral heavy mesons with mass E in HQET. Equation (17) shows that 
only depends on q 2 . Choosing a frame in which q = 0, we can use the analyticity of £„(g) 
to write a dispersion relation, 



v < ^ 1 [°° JF ImE,(£,0) 

Z v (m D ,0) = - dE— — . 

7i hm-n E — m D + ze 



(18) 



Using Eq. (17), we obtain 



Am = -—P dE K J +0 tlS£R) . (19) 
2vr J2m n [E-m D V E J\ y ' 

Eq. (19) is the main result of this section. It expresses Amo in terms of a weighted 
integral of the width difference of heavy mesons, AF(E), over varying heavy meson masses, 



2 In the OPE-based calculations, because m c /AQCD is not very large and sublcading terms in the Aqcd /m c 
expansion are enhanced by A X 5s/m s [10], such terms dominate the short distance contribution [10-12]. 
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E. The heavy quark limit was essential in deriving this relation, since has a physical 

interpretation for arbitrary q, while for q ^ p D , does not. The C(Aq CD / E) error 

in the integrand is a consequence of our reliance on this limit, and the resulting correction 
is 0(1) in the small E region. Dispersion relations for Am^ were considered previously in 
Ref. [13], where ImS(s) (with a different definition of S) was modeled, but it does not have 
a physical interpretation for s ^ m 2 D . 

To calculate x/y using the dispersion relation, we need to know Ar as a function of the 
heavy meson mass. Examining Eq. (19), we expect that values of E close to m D give the 
largest contribution to x. In the next section we recall the calculation of AT(E) performed 
in Ref. [3]. If AT(E) is a decreasing function of E at least as a positive power, 1/E a with 
a > 0, then the dispersion relation does not require subtraction in order to converge. In the 
model we consider, AT(E) actually falls off as ~ l/E 2 , and we will argue that some kind of 
decreasing behavior is likely to hold model independently. 

III. CALCULATION OF THE LIFETIME DIFFERENCE 

The computation of x using Eq. (19), requires us to know Ar for a heavy meson of 
varying mass. The calculation of Ar cannot at present be done from first principles. In 
Ref. [3] Ar was computed using a simple model in which SU(3) breaking was taken into 
account in calculable phase space differences, but neglected in the incalculable hadronic 
matrix elements. This approach was motivated by the fact that phase space differences 
alone can explain the experimental data in several cases; for example the ratio r(D| — > 
Dn)/T(D* 2 -> D*n) [14], the large SU(3) breaking in T(D -> K*lV)/Y{D -> plv) [15], and 
the lifetime ratio T£>Jt d o [16]. It certainly cannot explain all SU(3) violation, for example, 
T(D — > 7i7i)/T(D — > KK). The generic conclusion of Ref. [3] was that if multi-body final 
states close to the D threshold have significant branching ratios, then they can give rise to 
sizable contributions to Ar that are absent in the OFE-based calculations. Our purpose in 
the next section will be to see whether the same mechanism can also give rise to x at or 
near the percent level. Here we review the analysis of Ref. [3] . 

We denote a set of final states F belonging to a certain representation R of SU(3) by 
Fr. For example, for two pseudoscalar mesons in the octet, the possible representations for 
F = PP are R = 8 and 27. In Ref. [3] it was shown that uf r , the value which y would take 
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if elements of F R were the only channels open for D° decay, can be expressed as 





7~~C W 


n) (n\H w 


D°) _ EneF R (D° 




n)(n\H w 


D°) 


EneF R (D°\H w 


n) (n 


H W \D») En6Ffl r(L>o^ 


n) 





(20) 



The derivation of this relation assumes the absence of CP violation, so that (D°\Ti w \n) is 
related to {D°\H w \n), and uses the fact that both \n) and \n) belong to the same S77(3) 
multiplet. When the S77(3) breaking in the matrix elements is neglected, Eq. (20) gives a 
calculable contribution to yp R without any hadronic parameters. The numerator contains a 
combination of Clebsch-Gordan and CKM coefficients that ensures that df r is proportional 
to m 2 sin 2 6c when the sum over all members of any given multiplet Fr is performed, as 
required by Eq. (3). 

As an example, the contribution of the multiplet containing two pseudoscalar mesons in 
an SU (3) octet is given by 



V{PP)i 



= sin 2 0, 



c 



+ 1 $(77, tt ) - 1 $(77, K°) - 2$(K + , 7T- 



$(ir u ,7r' 



x 



1 , , 



1 _™ 



- $(?7, K°) + 7T+) + - $(K°, TT ) 



-1 



+ G(sin 4 ^), (21) 



where $(n) is the phase space factor for D — > n decay. Then y can be computed as the sum 
of the 2/f k 's weighted with the D° decay rate to each representation, 

!/=f E*4 E r(D°-n)j. (22) 

The were computed for all PP, PV, and W representations, and for the fully symmetric 
3P and AP final states [3]. The contribution of poles corresponding to nearby K resonances 
was shown to be small [3, 17]. Assuming that the values of y(4P) R for R = 8, 27, 27' are 
typical for all R, it was found that the 4P final states give a contribution to Ar at the 
percent level. The result is large because many of the decays in question are close to or 
above threshold, so the SU (3) cancellation in these multiplets is largely ineffective, yielding 
U(4P) R = 0(0.1) [3]. Moreover, the D° branching ratio to four pseudoscalars is approximately 
10%. 

We shall now use this model of 577(3) breaking, together with some assumptions about 
the energy dependences of the relevant decay rates, to compute x/y. 
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IV. CALCULATION OF THE MASS DIFFERENCE 



The crucial difference between the calculation of x and y is that once we assume that the 
only source of SU(3) breaking is from the final state phase space differences, the hadronic 
matrix elements cancel in y, but not in x. As determined by Eq. (19), x depends on AT(E), 
and so the ^-dependence of the hadronic matrix elements does affect x. Using Eq. (19), we 
find for x/y, 



_ xjn _ 1 p f 00 dE y FR (E) T Fr (E) ^ 
Vf r tt J2m 7r E -m D y FR (m D ) T FR (m D ) ' 



We will quote our results in terms of t Fr . To proceed further we need to understand or make 
some assumptions about the independence of the decay rate to the final state F, F F (E). 
We define the dimensionless function 

F$&- (24) 

and we will study the ^-dependence of this quantity. Note that the constant of proportion- 
ality in Eq. (24) cancels in the ratio rp R . Moreover, g F is expected to depend only on the 
final state F, and not on the SU (3) representation R. 

One can reconstruct x from x Fr using a relation analogous to Eq. (22). Below we calculate 
t Fr for several final states and then estimate the total x. First we will study F = PP, 
because it is a simple case that is interesting to understand in detail. Then we will turn to 
F = 4P, because it is the final state that can give y ~ 1%. 



A. Two-body D — ► PP decays 

For decays to two pseudoscalar mesons, it is possible to develop a reasonable model 
°f 9pp{E). When m H ^> A QCD , we may approximate the H — > -nit amplitude with its 
factorized form. Here A(H — > nn) ~ G f Vckm m 2 H f n Fn^n, where /„- is the pion decay 
constant and Fh^-k is the H — > n form factor at q 2 = m 2 . It has been shown that, as 
m H oo, Fh-tk oc (A/m H ) 3 / 2+x [18], where X arises from summing Sudakov logarithms of 
the form exp [Ca s {m H )\n 2 (m H /Kj\ ~ (A/m H ) x with X = -2nC/(3 . Since T oc \A\ 2 /m H , 
we conclude that 

g PP (E > Aqcd) oc E~ 2X . (25) 
The existing calculations suggest that \X\ <^ 1 [19], so we set X = hereafter. 
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FIG. 2: Predictions for r^ PP ^ s (solid curve) and r^ PP ^ 2T (dashed curve) as functions of mi. 

In the E — > limit our calculation is necessarily unreliable, as the derivation of Eq. (19) 
relied on HQET. Nevertheless, as a model we will take the behavior of the K — > tctt amplitude 
in chiral perturbation theory. At leading order, this transition is mediated by an operator 
of the form Tr^St (9<9 M £), where £ = exp[2iM / /] and M. is the meson octet. Since this 
term has two derivatives, it implies that the decay amplitude is proportional to m 2 K . Since 
this is the only dependence on ttik in the amplitude, the ^-dependence of the rate is 



where m 12 are free parameters. The overall normalization cancels in the results. This model 
allows for a "chiral" region, E < mi, an "intermediate" region, mi < E < m.2, and a "high 
energy" region, E > m<i- In our calculations mi is allowed to vary in the range 0.2 — 1.0 GeV, 
and 77i2 m the range 1.5 — 10 GeV. As we emphasized above, our derivation relies on HQET, 
so any strong dependence on scales below ~ 1 GeV would signal an irreducible lack of 
reliability. 

In Fig. 2 we plot ?"(pp) 8 (solid) and rr P p\ 27 (dashed) as a function of m 1? for m 2 = 2 GeV. 
In this case all members of the final state representations are kinematically allowed and have 
large phase space, so we find that the result is dominated by cancellations below the scale 
m D . Therefore r PP is sensitive to the shape of g PP (E) at low energies, i.e., the value of mi, 



g PP {E -> 0) oc E\ 



(26) 




(27) 
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but changing m 2 to 3 or 4GeV has little effect on r PP . Because of the strong dependence 
on mi, we should not trust this result. However, since y for these representations is very 
small, y(pp) R = —0.018% and y(pp) 27 = —0.0034% [3], these final states do not give sizable 
contributions to x in any case. 

When we consider decays to the lightest pseudoscalar octet, the dependence of these 
pseudo-Goldstone boson masses on m s is given by (for m u ^ = 0) 

m 2 n = 0, m 2 K = /im s , m ^ = \ ^ ms ' ^ 
where // is a hadronic scale. We can then expand AT(E) for large E as 



ATpp(E) = 



Tpp(E) 



Because SU (3) breaking in our approach comes from phase space differences, the coefficients 
Cj depend quadratically on the masses of the final state particles. Since in Eq. (28) m s is 
always accompanied by \i and Ar must be suppressed by m 2 s , we conclude that Co = c\ = 0. 
The coefficient c 2 can be proportional to /j, 2 m 2 and is the leading nonvanishing term, implying 
a l/E 4 suppression of ATpp(E) compared to T PP (E). However, the actual it, K, and r\ 
masses do not exactly satisfy Eq. (28) in the m U:d = limit, nor the Gell-Mann-Okubo 
(GMO) relation, 3m 2 = Am 2 K — m 2 . Violating the GMO relation is equivalent to adding 
a small term to m 2 K or to m 2 of the form em 2 . This changes the asymptotic behavior 
of AT(E), because now we can have C\ ~ em 2 . Since the D — > PP decay is far from 
threshold, the 577(3) cancellation in this channel is very sensitive to the pseudoscalar meson 
masses. This can be verified analytically by expanding Eq. (21). As shown in Fig. 3 (again 
for m 2 = 2GeV), imposing the GMO relation on the n, K, and r\ masses decreases r PP 
significantly, in such a manner that ypp increases by roughly the same factor, while \x P p\ 
is approximately stable at the (5 — 8) x 10~ 4 level. As discussed in Ref. [3], our results have 
little sensitivity to including or neglecting 77 — 77 — 77' mixing. 

By contrast, for final states including vector mesons or heavier pseudoscalar represen- 
tations, the masses of the mesons depend linearly on m s . Thus, for these final states, 
AT F (E) /V F (E) is simply proportional to m 2 /E 2 for large E, and there is no strong de- 
pendence on the precise values of the hadron masses. This is the minimal suppression of 
AT F (E) /V F {E) consistent with group theory, i.e., Eq. (3), and our phase space model for 
S77(3) violation indeed gives such an effect. These results imply that the dispersion relation 
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FIG. 3: Predictions for ?"(pp) g (solid curve) and ?7pp) 27 (dashed curve) as functions of mi, imposing 
the GMO relation on the tt, K, and r\ masses. 

in Eq. (19) converges for any final state F, for which Ff(E) does not increase as E 2 or 
faster. This is very likely to be true for all final states (recall that Tpp(E) ~ constant for 
large E). 



B. Four-body D — > 4P decays 

Now we turn to the 4P final state in the fully symmetric 8, 27, and 27' representations 
of SU(3). We know even less about g±p(E) than about g PP (E), so we use two models to 
attempt to bracket roughly the uncertainties, 

{Ejm\ for E < vii\ , 
1 for m 1 < E < m 2 , (30) 

m2 /E for E > m 2 . 

The choice of g' 4P (E) allows for the possibility that T(H — > 4P) may start to fall for large 
ran instead of remain constant. This alternative is motivated by the argument that because 
the quasi-two-body picture holds only in a small part of phase space, in most of the phase 
space the opening of many decay channels will reduce the rate. 

The left plot in Fig. 4 shows r( 4 p) g (solid curve), T( 4 p) 27 (long dashed curve), and ?"(4P) 27 , 
(short dashed curve), as functions of m 2 , using g±p(E) with mi = 0.8 GeV. For m\ <1 GeV 
there is no dependence on m\. The dependence of the curves on 777,2 is negligible for 7772 ^ 
3 GeV. If we use g' 4P (E) instead, shown in the right plot in Fig. 4, then r^ P ) R changes 
roughly by a factor of two. We have explored other forms of g±p(E) as well, and we find 
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FIG. 4: Predictions for r( 4 p) g (solid curve), 7"( 4 p) 27 (long dashed curve), and r(4P) (short dashed 
curve), as functions of m2 for the models gip{E) (left figure) and g' iP {E) (right figure) in Eq. (30). 

that these two cases cover a reasonable range of predictions. 

In contrast to D — ► PP decays, for the 4P final state there is no strong dependence on 
the 7T, if and 77 masses. Because the decay is close to threshold, the dispersion integral is 
dominated by E near mo, where some of the 4P final states are kinematically forbidden, 
and so the sensitivity to the pseudoscalar meson masses is reduced. Imposing the GMO 
relation makes only a small difference; for example, for the (4P)§ representation the value 
r( 4 p) g = —0.98 obtained with the g± P (E) model, mi = 0.8 GeV, m 2 = 3 GeV, and the 
physical meson masses (corresponding to the solid curve in the left plot in Fig. 4), would 
change to r( 4 p) g = —0.87 if the GMO relation were imposed. 

V. DISCUSSION AND CONCLUSIONS 

It is likely that the dominant contributions to the mass and width differences in the D 
system have a long distance origin in the SM. Therefore, naively one would expect x and y 
to be of the same order of magnitude. We have derived a new dispersion relation (19) and 
used it to study this question. Our dispersion relation has the useful property that it relates 
the mass difference in the heavy neutral meson system at fixed heavy meson mass to the 
physical width difference of heavy mesons with varying mass. 

The advantage of using a dispersion relation that relates x to y is that we can use 
existing models for y to calculate x. Our dispersion relation is likely to converge without 
any subtraction, because the SU(3) breaking required to yield nonzero mixing introduces 
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an m 2 s /E 2 suppression in y(E). We have used a model in which SU(3) breaking arises from 
phase space differences, which may give a reasonable approximation to y(E) only when E 
is not very large. Since the derivation of the dispersion relation employed the heavy quark 
limit, it is essential not to interpret our analysis as a precise calculation for x. Instead, we 
used this model only to get a rough and qualitative prediction about the likely relation of x 
to y. 

To make numerical predictions we needed the heavy mass dependence of heavy meson 
partial widths to certain final states, which introduces some additional model dependence 
in our results. (For decays to two pseudoscalars, there are limits in which one can draw 
firmer conclusions about the mass dependence, which we have incorporated into the model.) 
We calculated the ratio x/y for PP and 4P final states. Our conclusion is that it is indeed 
likely that in the Standard Model, x is not much smaller than y in the D system. In our 
numerical study, we found that for the 4P final state, x/y varies roughly between —0.1 and 
— 1. We conclude that if y is in the ballpark of +1% as expected if the 4P final states 
dominate y [3], then we should expect \x\ between 10~ 3 and 10 -2 , and that x and y are of 
opposite sign. This estimate has a large uncertainty, and we can trust it only at the order 
of magnitude level. We have explored the sensitivity of this qualitative result to a number 
of the assumptions we have made, and have found that changing the details of the model 
does not significantly alter our conclusions. Furthermore, including some SU (3) breaking in 
the matrix elements cancels to some extent in x/y and does not induce dramatic changes. 

The significance of our result is clear only in the context of the experimental situation. 
The current bounds on x and y are at the level of a few percent, and the central question is 
whether their actual observation at or just below this level could be interpreted as a clear 
signal of physics beyond the Standard Model. We would argue that our analysis has taught 
us that, without further refinement, the answer is no. We have identified a real effect that 
could plausibly give x and y at the percent level, albeit with very large uncertainties. 

In general, an observation of x ^> y would be an indication for new physics, but this 
could only be established if y were very small, at the 10~ 3 level. Such a situation could 
arise if new physics enhanced x but not y. Yet since one cannot exclude the possibility of 
cancellation between different SM contributions to y, even this outcome would not admit 
an unambiguous interpretation. 

However, if x were indeed enhanced by new physics, such new physics may also introduce 



14 



a sizable new CP violating phase which may be observable. Thus, we would argue that in 
D° — D° mixing, the only single measurement that could establish by itself the presence of 
new physics would be the observation of CP violation, which is very small in the Standard 
Model independent of hadronic effects. 
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